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We present a new method to measure FOj in flux-limited samples of active galactic nuclei (AGN) 
jets, where F is the bulk Lorentz factor and Oj is the jet's half-opening angle. The FOj parameter 
is physically important for models of jet launching, and also determines the effectiveness of jet 
instabilities and magnetic reconnection. We measure FOj by analyzing the observed distribution 
of apparent opening angles in very long baseline interferometry (VLBI) flux-limited samples of 
jets, given some prior knowledge of the active galactic nuclei (AGN) radio luminosity function. 
t-H We then apply this method to the MOJAVE flux-limited sample of radio loud objects and find 

FOj « 0. 1 ± 0.03, which implies that AGN jets are subject to a variety of physical processes that 
require causal connection. 
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1. Introduction 

A relativistic jet's value of the parameter Fdj can impact jet physics for several different rea- 
sons, where T is the bulk Lorentz factor and dj is the jet's half-opening angle. In some jet models, 
Fdj is expected to be ~ 1 if either the jet's Poynting flux is efficiently converted to kinetic flux [1], 
or if the jet is freely expanding. For hydrodynamic models, Fdj is related to the ratio of the pressure 
of the ambient medium to the jet's internal pressure [2]. Also, most jet acceleration models require 
that the jet be in causal contact, which implies that Fdj < 1. Besides global jet dynamics, Fdj may 
also have significance for transient events such as magnetic reconnection and jet instabilities, since 
the jet must be in causal contact for such processes to occur [3], and Fdj in part controls how much 
jet sidewise expansion slows instability growth [4, 5] . Additionally, Fdj is a significant parameter 
affecting the structure of shocked boundary layers [6] and the efficiency of particle acceleration at 
recollimation shocks [7] . 

There have been two past measurements of the characteristic value of Fdj for AGN jets. Using 
7mm very long baseline array (VLBA) data from 15 different AGN jets, Jorstad et al. [8] measured 
Fdj by assuming that the observed pattern speed of moving jet components corresponds to the jet 
bulk flow speed, and that the component variability times are equal to the jet frame light crossing 
times of the (resolved) components. From these assumptions they determined the component's 
Lorentz factor and viewing angle, and found an anti-correlation between the derived values of dj 
and T, with Fdj = 0.17 ± 0.08. With a larger sample of 56 AGN jets from 15 GHz VLBA data, 
Pushkarev et al. [9] performed the same analysis, except they used Lorentz factors and viewing 
angles from Hovatta et al. [10], who derived these values used variability time, maximum flux 
density of flares, and equipartition derived brightness temperature arguments. Pushkarev et al.'s [9] 
analysis found a similar anti-correlation with Fdj = 0.13. 

Motivated by the above theoretical concerns, we set out here to provide a new method of 
measuring Fdj by using a larger sample of AGN jets and very different physical assumptions. We 
do this by deriving the expected distribution of apparent opening angles in a flux limited sample 
with Fdj as a free parameter to be fixed by finding the best fit to an empirical distribution of 
apparent opening angles. Our data consist of the 135 AGN jets that make up the MOJAVE-I sample, 
a 15 GHz flux density-limited survey conducted by the VLBA of radio sources in the northern sky 
with flux densities above 1.5 Jy, and above 2 Jy for sources with —20 < dec < [11]. 

2. Constraining AGN jet physics with observed apparent opening angles 

Here we derive the probability density function (PDF) for apparent opening angles, P(0 app ), 
and then fit it to a distribution of apparent half-opening angles, app , measured from MOJAVE-I 
data. To test the viability of the simplest case scenario, we assume that Fdj is constant for all 
relativistic AGN jets. As we show below, Fdj is then a free parameter of P(d app ), and thus will be 
determined in the fit to the empirical distribution of app . 

To derive P(d app ), we first derive the PDF for viewing angles, P(d). If the sample of AGN jets 
were unbiased with respect to orientation, P(d) = sin d. However, because it is flux-limited, it will 
take the form 

P(d) = Doppler bias factor x sin0. (2.1) 
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This additional factor takes into account that more sources are directed at the observer in a flux 
limited sample because Doppler beamed jets are detectable at further distances than unbeamed 
ones. Vermeulen & Cohen [12] computed this term, finding that it depends on the bulk Lorentz 
factor distribution in a flux limited sample (every jet is assumed to possess a single Lorentz factor), 
the integral source count index, and the beaming index of the jet. The beaming index is denned 
from the relation F = 8"F', where n is the beaming index, 8 is the Doppler factor, F is the observed 
flux density, and F' is the intrinsic flux density. The observed integral source count index is defined 
in the expression N(> F) oc F~ q , representing the number of sources N observed with a flux density 
above F, which is a power law in F with source count index q. Including the Doppler bias factor 
in the viewing angle PDF gives 

P(6,F) =A(l- J 8cos0)^ 1 sin0P(r), (2.2) 

where a = nq — \,P(T) is the PDF for jet bulk Lorentz factor, and A is the constant of normalization. 
An important assumption made in calculating the Doppler bias term is that the log-log slope of 
N{> F') vs F' and N(> F) vs. F are the same, which Vermeulen & Cohen [12] justify based on 
previous studies of AGN jet luminosity functions. 

The opening angle distribution may now be derived from P(6,T) by a change of variables 
from 6 to app , and marginalizing over T: 



p(e app ) = J drp(e(d app ,r),r) 



do 



(2.3) 



where the 6 and dd/dd app are functions of app and T, and can be determined by assuming a 
particular jet geometry, which we take to be conical here. These relationships are often derived 
by treating conical jets as triangles projected onto the plane of the sky, implying that tan app = 
Rj/C = Rj/ (£sin0), where 8 is the jet viewing angle, Rj is the jet radius, I is the jet length, and £' 
is the jet length projected onto the sky. If we assume dj Rj/£, then 

tan dj = tan app sin 6 . (2.4) 

For app , dj <C 1, this reduces to the most commonly used relation for jets, dj = app sin0 [8, 9]. 
However, this relation breaks down as 6 — > dj, since when 6 = dj, equation (2.4) incorrectly gives 
$app = ft/4 (it should be 7r/2 in this case because the projected jet streamlines would cover half of 
all directions on the sky). However, we ignore such situations here because most jets are viewed 
such that 6 > 6j [13]. 

We now calculate an approximate analytic expression for P(0 app ) and show that it is not very 
sensitive to P(T). Note that in a flux-limited VLBI sample, jets with small viewing angles will 
dominate, so we assume 0<1 and F> 1, and approximate (2.2) as 

p ( ' r )=^^pr(i+r 2 e 2 )- fl - 1 e. (2.5) 

For simplicity, we now evaluate equation (2.3) in light of the geometry implied by equation (2.4), 
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Figure 1: Fitting apparent opening angle data to the derived probability density (eqn. 2.6) via maximum 
likelihood estimation. The only free parameter is p, and a = 3 is assumed. The result is p fit = 0.095 ± 0.009. 
Note that the best fit cumulative distribution function (CDF) is systematically lower than the empirical CDF 
for 6> app ;> 0.2. 

with the additional approximation that sin0 « 6, and obtain 



where p = F8j, and we have absorbed the term in brackets into the new normalization, A'. Thus, 
it is apparent from equation (2.6) that P(d app ) does not depend significantly on the form of P(T). 
Equation (2.6) is an accurate approximation of P(6 app ) as long as p <C 1, which, as shown below, 
is a valid assumption. 

We now fit equation (2.6) to the empirical opening angle distribution from the MOJAVE-I data 
set using maximum likelihood estimation (MLE), as shown in figure 1. We set a = 3, so our fits 
only include one free parameter, p (= T0j). This value of a is close to the value it would have if the 
beaming index corresponded to a steady jet [14] with a spectral index of a = 0.7 (F v °< v~ a ) [15], 
and q = 1.5 [12]. The result we obtain is pfi t = 0.095 ±0.009 using equation (2.4). The confidence 
interval, calculated from the inverse of the information matrix, is not very meaningful given the 
uncertainty in a discussed below and other theoretical uncertainties discussed in §3. 

The true confidence interval in the best fit value of p requires analyzing a, whose theoretically 
likely values range from ~ 2 to 5. This range originates from beaming models of jets, which bound 
the beaming index n as being between 2 + a and 3 ± a, and the measured integral source count 




(2.6) 
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index of q « 1.5 [12]. The spectral index, a, for jets is typically between and ~ 1 [15], thus 
given the range of n as 2 + to 3 + 1 , we find the range of a as being between 2x1.5 — 1= 2 and 
4 x 1.5 — 1 = 5. Assuming such a range in a alone determines the confidence interval of the best 
fit, then p fit = 0.096 ± 0.025. 

To better understand how best fit values of p depend on the assumed value of a, we calculate 
the best fit value of p for different given values of a. We find that the resulting curve of Pfi t versus 
a is very well reproduced by assuming that each best fit value of p is sensitive to the location of 
the peak of the empirical number density distribution, 0a Pp k ~ 0.1 to 0.2, which is the inflection 
point of the cumulative distribution shown in figure 1 . This assumption implies that the relationship 
between p and a is can be found by maximizing the probability density P(0 app ), which according 
to equation (2.6) gives p m 0jpp \J1a—\, or 

a peak 

TO, re - " app V - (2.7) 

1 y/2{nq-l)-l 

For an empirical peak of 0fp P = 0.13, this correctly reproduces the best fit values of p for values 
of a between 2 and 10 to within an error of 7%. 

Notably, our best fit theoretical cumulative distribution function (CDF) does not reproduce the 
empirical CDF for app ^ 0.2. This is not surprising giving the simplicity of our one parameter 
model, and that large apparent opening angle sources are more likely to have smaller viewing 
angles 6. Sources with small 6 are more subject to differential Doppler beaming, where different 
streamlines in the jet have significantly different Doppler factors. Other effects our simple model 
does not take into account, and that might explain the deviation between our data and best fit, are 
discussed below in the conclusion. 



3. Conclusion 

We have illustrated a new method that measures the characteristic value of Tdj in a flux- 
limited sample of jets, yielding Tdj wO.l ±0.03. This measurement implies that AGN jets are 
causally connected and therefore subject to a variety of processes requiring causal connection, 
such as large-scale magnetic reconnection and the development of global instabilities including the 
current-driven kink mode, pressure-driven modes, and Kelvin-Helmholtz modes. 

Some caution is needed in interpreting these results due both to our model's simplicity and the 
deviation of our best fit theoretical CDF from the empirical CDF for large apparent opening angle 
sources. The derived PDF used to fit the empirical distribution of apparent opening angles crucially 
depends on the geometry of the emitting region, which we naively assumed to be a simple cone. 
This is likely not the case in jets subject to relativistic velocity shear, where differential Doppler 
beaming causes different parts of the jet to become observable depending on the jet viewing angle. 
Thus, the equation for how apparent opening angle depends on jet viewing angle and intrinsic 
opening angle (eqn. 2.4) may need to be modified to reflect the effect of velocity shear. Also, the 
assumptions implicit in Vermeulen & Cohen's [12] calculation of the viewing angle PDF are open 
to criticism. Namely, the result that viewing angle PDF only depends on the slope of the log-log 
integral source count plot, and not other parameters of the AGN jet luminosity function, needs to 
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be further investigated. We intend to address the above issues in a future publication. Despite the 
above issues, it is encouraging that our results of Fdj « 0. 1 ± 0.03 so closely aligns with previous 
measurements of Fdj = 0.13 [9] and Tdj = 0.17 ±0.08 [8], which used very different methods and 
assumptions. 
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